Accepted for publication in ApJ Letters, February 15, 2010 

Preprint typeset using BTpX style emulateapj v. 1 1/10/09 



DIRECT OBSERVATION OF A SHARP TRANSITION TO COHERENCE IN DENSE CORES 

Jaime E. Pineda 1 , Alyssa A. Goodman 1 , Hector G. Arce 2 , Paola Caselli 3 , Jonathan B. Foster 1 ' 4 , Philip C. Myers 1 , 

Erik W. Rosolowsky 5 

Accepted for publication in ApJ Letters, February 15, 2010 

ABSTRACT 

We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). 
The map covers a region large enough (~ll'x 14') that it contains the entire dense core observed in previous 
dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum 
features found in bolometer observations. The velocity dispersion in the central region of the core is small, 
presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage 
and high sensitivity of the observations that we present the detection, for the first time, of the transition between 
the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the 
velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). 
The change in velocity dispersion at the transition is « 3 krns" 1 pc -1 . The existence of the transition provides 
a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. 
From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition. 

Subject headings: ISM: clouds — stars: formation — ISM: molecules — ISM: individual (Perseus Molecular 
Complex, B5) 



1. INTRODUCTION 

The velocity dispersion in molecular clouds (MCs) has been 
known for years to be supersonic. Several numerical simu- 
lations of supersonic turbulence can successfully reproduce 
some of the MCs properties. But, dense cores, where stars 
are actually formed, present velocity dispersions with non- 
thermal motions small er than the thermal values and also 
independent of scale (Goodman et al. |1998 Caselli et al. 



2002). |Goodman et aL] ( |1998| ) anaiTaselli et al.|( |2002| ) coined 
the term "coherent core," to describe the region where the 
non-thermal motions are subsonic and constant as " i slands 
of calm in a more turbulent sea." [Goodma n et al.| ([1998) 
showed that the lower density gas around cores, traced by 
OH and C 18 (1-0), presents supersonic velocity dispersions 
that decrease with size, as expected in a turbulent flow, while 
the dense gas associated with cores traced by NH 3 shows a 
nearly-constant, nearly-thermal width. Therefore, a transition 
between turbulent gas and more quiescent gas must happen at 
some point. However, it was not clear if the transition could 
be detected in the dense gas tracer on its own and/or if it is a 
smooth or abrupt transition. 

The nearby Perseus MC, at ~250 pc, is a good place to 
search for the transition to coherence. A large number of 
' 'dens e cores" have been identified in dust continuum sur- 
veys ( [Hatchell et al.|[2005] |Enoch et al.| [20061 |Kirk et al. 



2006), and previous NH 3 observations have been carried out 



Rosolowsky el^ ( [200% |Ladd et al.| ( [T994| ; |Benson & My 
ers|(|1989|). From these core lists, B 5 stands out as a rather 
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isolated an d bright dense core that has been studied in detail 
in the past (Young et al.|1982[|Langer et al.|19891|Fuller et al. 
199T| |Bally et al.||1996[ |Bensch||2006[ |Pineda et al.| |2008 ), 
with ~ 1 1 Mq detected in dust continuum. It should be noted 
that dense cores are usually identified as objects which are de- 
tected in some molecular line emission tracing dense gas, e.g., 
NH 3 (1,1) orN 2 H + (l-0), or in dust continuum, e.g., MAMBO, 
SCUBA, BOLOCAM, LABOCA. But, it is not clear that dif- 
ferent tracers and instruments provide a consistent object def- 
inition. 

If velocity dispersion and density profiles are related, then 
a transition should also be seen in column density. For this 
at least two well known techniques could be used: extinction 
maps and dust continuum emission maps. Extinction maps 
trace the total column density along the line of sight, an d they 
are useful to study the large scale structure in MCs (e.g . , |Lada 



[ et al.|1994t|Alves et al.| 1998 1 Lombar di^lves|2001t|Pmeaa 
|et al.|2008[ |Goo dman et al. 2009) and to identify dense cores 
(e.g., |Alves et al.||2007| ). Unfortunately, the coarse angular 
resolution achieved in extinction maps (> l r ) is still not quite 
enough to allow detailed studies of the column density struc- 
ture within and around c ores, with the ex ception of only a 
handful of cases e.g., B 68 ( [Alves et al. 2001 ]) Dust continuum 
obser vations usin g SCUBA (|Hatchell et al.||2005||Kirk et al.| 
2006), MAMBO (Motte et al.||1998[ Kauffma nn et aL|f 2008) 
and SHARC II ( |Li et al.|200'7] ) have angular resolution oTL5*, 
11" and 9" , respectively, which are high enough to enable the 
study of the transition to coherence. However, due to observ- 
ing techniques (sky removal and/or chopping), they filter out 
large-scale emission (usually between 1.5' and 2'), and there- 
fore dust continuum maps are not suitable to study the region 
where the transition to coherence happens. Currently, these 
limitations on dust mapping leaves high-resolution, spatially - 
unfiltered, molecular line observations as the best tool to look 
for the "environs" to "core" transition. 

In this letter, we present new NH 3 observations of B5 ob- 
tained with the 100-m Robert F. Byrd Green Bank Telescope 
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Fig. 1 . — Left panel: Map of 13 CO(1-0) integrated intensity for the B5 region obtained by the COMPLETE Survey ( Ridge et al. 2006}. The orang e box shows 
the ar ea mapped with the GBT. Right panel: Integrated intensity map of B5 in NH3(1,1). White contours are BOLOCAM dust continuum emission (EnochetaT] 



|2006) , while gray contours show the 0.15 and 0.3 Kkm 
direction is shown by the arrows. 

(GBTQ which provide the first detection of the transition to 
coherence in a single tracer. These observations provide an- 
swers to two questions: a) what is the extent of coherent dense 
cores?; and b) is the transition to coherence smooth or abrupt? 

2. DATA 

We observed B5 using the GBT. The observations were car- 
ried out between December 23 and March 31, 200 9 (project 
08C-088) , using the On-The-Fly (OTF) technique (Mangum 
et al. |2007] ), with a dump rate of 3 dumps per beam, and pro- 
ducing a dump every 3 seconds. 

We used the high-frequency K-band receiver and config- 
ured the spectrometer to observe four 12.5 MHz windows cen- 
tered on NHsCU), NH 3 (2,2), CCS(2i-l ) and HC 5 N(9-8) 
rest frequencies. We chose to use two feeds and two polar- 
izations simultaneously, which trades decreased spectral res- 
olution for increased sensitivity given GBT spectrometer con- 
straints. The spectrometer generated 4096 lags across each 
window, giving a 3.050 kHz channel separation, equivalent to 
0.04 krns" 1 for the NH 3 spectra. We observed in frequency 
switching mode, with a shift of 2.0599365 MHz around the 
center of the band. This configuration ensured that all 18 
hyperfine components of NH3(1,1) were observed within the 
spectral window. The pointing model was updated every 
60-90 minutes, depending on weather conditions, using the 
quasar 0336+3218. Flux calibration was carried out by ob- 
serving the flux calibrator 3C48 during each session. All 
the intensities reported here are on the 7^* scale, which is es- 
tablished using atmospheric opacity estimates at 22-23 GHz. 
Data cubes are generated using all observations taken and 
convolved onto a common grid with a tapered Bessel func- 
tion (see Mangum et al.|[2007] ). The GBT main beam effi 



level in NH3(1,1) integrated intensity. The young star, B5-IRS1, is shown by the star. The outflow 



ciency (rj m b) is 0.81 at these frequencies. All the data reduc- 
tion was car 
is 0.046 K. 



requ 

tion was carried out in GBTIDLQ The median rms in the map 



The resulting NH3(1,1) integrated intensity map for B5 is 
shown in Figure [T] and it covers a region of size ll'xW. 
Gray contours in Figure [T] show the extension of NH 3 (1,1) 
emission. 

3. RESULTS 



° F ull description and analysis of all of the COMPLETE Survey's (Ridge 
|etal.|2006) GBT NH 3 maps of Perseus cores, including B5, will be presented 
in Pineda et al., in preparation 
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Fig. 2. — Centroid velocity map of B5 obtained by fitting the NH3 lines. 
The position of the protostar is shown by the star. Black contours are Bolo- 
cam dust continuum emission (Enoch et al. 2006). The GBT beam size is 
shown at the bottom left. 

The NH3(1,1) and (2,2) li nes are fitted simu l taneou sly us- 
ing a forward model as in Rosolowsky et al.| ( |2008| ). This 
method allows us to obtain centroid velocity (vlsr), velocity 
dispersion (cr v ), kinetic temperature (71), excitation tempera- 
ture (T^) and opacity (rn) for every position, while also in- 
cluding the response of the frequency channel using a sine 
profile. If the NHs(l,l) line is optically thin then rn and T ex 
can not be obtained independently, and therefore the optically 
thin approximation is used if rn < 1 or a Tn < 0.4 rn, where 
a Tn is the optical depth uncertainty (obtained from the fit). 
The centroid velocity map is shown in Figure [2] for positions 
where NH 3 (1,1) is detected. When compareato dust emis- 
sion maps from BOLOCAM (see contours in Figures [T]and|2]) 
or SCUBA, we find that the NH3(1,1) emission is spatially 
more extended than its dust continuum counterpart. There- 
fore, dense gas traced by NH3(1,1) is detected outside the 
boundaries of the dust-defined "dense core," calling into ques- 
tion the accuracy of dense core classification based only on 
the detection of a high-density tracer. 

Since our subsequent analysis of a v depends on very high 
accuracy, we eliminate from further consideration positions 
that do not fulfill the following criteria: a) clear detection of 
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Fig. 3. — Velocity dispersion map derived from fitting all hyperfine components simultaneously. The protostar position is shown by the star, and the contour 
shows the contour T pea ^ = 0.5 K. The box on the map presents the region where we zoom in and present the corresponding spectra in the right panel, which shows 
only the main component of the NH3(1,1) line. The centroid velocity and velocity dispersion obtained from the fit is displayed for each position. Top spectra 
in the right panel display two main hyperfine components clearly separated thanks to their low velocity dispersion (the coherent core), while when moving to 
positions outside the core the lines get weaker and broader (evident by the disappearance of the gap between hyperfine components). 



NH 3 (1,1); b) a Gv < 0.05 kms" 1 ; and c) a Gv < 0.2<r v ; where 
a Gv is the a v uncertainty from the fit. These conservative 
criteria eliminate observations where the velocity dispersion 
are poorly determined. The velocity dispersion map for B5 is 
presented in Figure [3] From this map it is clear that the cen- 
tral part of the core presents a region with small and uniform 
velocity dispersion. However, outside of this region there is 
extended NH3(1,1) emission with much larger velocity dis- 
persion, as large as a four times the velocity dispersion found 
within the core. Therefore, we have observed/<9r the first time 
in a single tracer the transition between dense but turbulent 
gas into more quiescent dense gas. A more detailed view of 
the transition is shown in the right panel of Figure [3] which 
displays only the main hyperfine blend of NH3(1,1) along a 
vertical cut marked by the vertical box in the left panel, and 
where the gap between hyperfine components disappear for 
broader lines. The transition appears to be sharp at GBT res- 
olution, occurring over just one beam width (two 15.5" pix- 
els). Using a 250 pc distance to Perseus, the 31" beam gives 
a 0.04 pc upper limit to the transition scale. 

To better characterize the transition between turbulent and 
calm gas, we map out the local gradient in cr v over B5. The 
distribution of this gradient's absolute value, |Va v |, is shown 
as an inset histogram in Figure [4] (Note that only positions 
marked by arrows in the Figure's grey scale map of a v are in- 
cluded in the histogram.) The orientation of Va v at the tran- 
sition is almost perpendicular to it, while inside the coherent 
core it is close to randomly oriented. The gradient ampli- 
tude distribution exhibits two peaks, with a narrow peak at 
small dispersion coming from within the coherent core. For 
the region in and immediately surrounding the coherent core 
in B5 (where arrows are shown in Figure]?]), a typical value 
of 3 kms _1 pc _1 for |V<r v | is found. Combining this typical 
value of | Va v | with the difference between the median veloc- 
ity dispersion for positions with subsonic and supersonic non- 
thermal motions (0.32 kms -1 and 0.13 kms -1 , respectively) 
we estimate a transition physical scale of 0.06 pc. This es- 



timate of the transition scale is actually also an upper limit, 
because it does not take into account that the observations are 
smoothed by the telescope beam (0.04 pc at the distance of 
Perseus). 

In the right panel of Figure]?] we show the velocity disper- 
sion as a function of peak antenna temperature (T pea k). Points 
marked in red are at a distance smaller than 63" from the em- 
bedded YSO in B5 and likely have increased velocity disper- 
sion as a result. If the peak antenna temperature is used as 
a proxy for the distance from the core ce nter (as in |Barranco 
|& Goodman|1998[ |Goodman et al.|1998| ), then it is clear that 
close to the center of the core velocity dispersions are small 
and display a small spread (the coherent zone). At lower T pea k 
(larger radii) there is a sudden increase in a v . Notice that 
the uncertainty in the dispersions is comparable to the sym- 
bol size at T pea k > 0.7 K, and still relatively small even at 
the lowest intensities analyzed here. To re-assure ourselves 
that there is no bias in our fitting toward finding higher dis- 
persion for weak lines, we performed tests on synthetic data, 
and found no bias that could explain the trend in Figure [4] 
In fact, because the integrated intensity (ex T pea k^v) map is 
smooth (most likely due to a smooth column density profile) 
T pea k must rapidly decrease to compensate for the sharp tran- 
sition in a v . This very simple argument can explain the effect 
seen in Figure [4] however, it does not provide an answer to 
the origin of the velocity dispersion transition. 

4. DISCUSSION 

The detection of a sharp transition to coherence provides 
very stringent constraints on numerical models of dense cores. 
Certainly the study of the density structure is important to un- 
derstand the relation between the core and its environment, 
and also to stu dy the relation be tween density and velocity 
dispersion (e.g. |Myers & Fu ller 1992), however, such a dis- 
cussion is beyond the scope of this letter. Here, we present a 
transition in velocity dispersion, for which we can not confirm 
nor rule-out an analogous density transition. 
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Fig. 4. — Left panel: Velocity dispersion map in the background. The young star, B5-IRS1, is shown by the star. The velocity dispersion gradient is shown 
by the arrows. The inset shows the velocity dispersion gradient distribution. The narrow peak at low values are the positions within the coherent core, while the 
transition surface populates the second distribution of points at larger values. Right panel: Velocity dispersion as a function of peak antenna temperature (J peak)- 
Red points are positions with a distance to B5-IRS1 smaller than 2.025 times the GBT beam. The peak antenna temperature can be used as a proxy of distance 
from the center, where high T pea k are found close to the center of the core and low T pea ^ are found are larger distances. Notice the small dispersion in cr v at high 
Tpeak (the coherent core) consistent with subsonic non-thermal components, and the abrupt increase in a v when approaching the transition to coherence with 
supersonic non-thermal components. The dispersion in cr v within the coherent core is even smaller after considering that the red points are affected by the central 
YSO. Blue horizontal lines show the expected velocity dispersion for two values of the velocity dispersion non-thermal component ((tnt)'- 0.5c S5 ave, and c s ,ave, 
where c s ,ave is the sound speed of the average particle (fi = 2.33) assuming Tkj n =10 K. 



The presence of the sharp transition allows for a robust def- 
inition of a coherent dense core: a region with nearly-constant 
subsonic non-thermal motions. Most certainly, the proposed 
approach of using the transition to coherence to define a dense 
core is not as time-efficient as using only large format bolome- 
ters, but it provides an identification system that is based on 
a physical quantity, and therefore it should be consistent with 
more sensitive observations. In the future, when more obser- 
vations of the transition to coherence in molecular lines are 
available for cores also mapped in dust, it might (or might 
not!) be possible to develop an empirical relation to improve 
the coherent cores identification using only dust maps. 

The velocity dispersion cumulative distribution is shown in 
Figure [5] The transition to coherence is a distinct feature in 
the cumulative distribution, where a change in slope is clearly 
observed. The effect is not only local, but it is also evident 
in the cumulative distribution for the entire region. Moreover, 
the velocity dispersion at which the cumulative distribution 
slope changes is robust against variations in the region se- 
lected to generate the cumulative function. 

A study comparing kinetic temperature and velocity disper- 
sion across the transition would be important to understand its 
origin. However, outside the coherent core, increases in ve- 
locity dispersion are accompanied by decreases in line bright- 
ness. In our present GBT data set, the NH3(2,2) emission 
beyond the transition to coherence cannot be reliably mapped 
because of the weaker lines. Since both (1,1) and (2,2) mea- 
surements are needed to determine kinetic temperature, we 
c an not yet study how temperature varies across the transition. 

lAlves et aL] ( |2001| ) showed that the column density profile 
of B68 can be well modeled by a Bonnert-Ebert (BE) sphere. 
Since then, this analysis has been applied to more cores find- 




FiG. 5. — Velocity dispersion cumulative distribution. Red curve shows the 
cumulative distribution for all points with good velocity dispersion measure- 
ments, while the black curve uses only points close to the central core. The 
sharp transition in velocity dispersion produces a change in the cumulative 
distribution's slope, which can be observed both locally (central region) and 
globally (entire map). The inset shows the spatial distribution of the points 
used, and velocity dispersion contours are overlaid in gray. 



ing that it usually is a good description (e.g., Ka ndori et al. 



2005), although we do not try to model B5 as a BE sphere. 
However, a column density profile simila r to BE can also be 
obtained in more dynamic eve nts (e.g. |Myers||2005} |G6mez 
et al. 12007] ). |Lada et al.| ( |2008| ) argued that most of the cores 
in the P ipe can be pressure confin ed b y the MC's own weight , 
see also |Bertoldi & McKee| ( |l992| ) and |Johnstone et a l. ( 2004). 
The observed increase in the velocity dispersion might be evi- 
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dence for a pressure difference between the coherent core and 
the external medium. However, in all these cases there is no 
explanation or description of what happens at the core bound- 
ary: is there a discontinuity? or is it a smooth transition with 
the background? 

The presence of a sharp transition to coherence suggests 
shock and/or instability /fragmentation origins. Shocks are 
predicted in models of core formation in supersonic flows 



(Padoan et 



m et al.||1997|) and in mod els of colliding large-scale 
e.g. |Heitsch et alT] ([2005 ). Core formation simula- 



nows, 

tions (ID) from co nverging supersonic flows ( [Gomez et al. 
|2007[ |Gong & Qstr iker 2009 ) predict a density and velocity 
discontinuity at the (isothermal) shock front position, which 
would also provide a core definition. Unfortunately, there is 
no discussion of th e spatial dependenc e of the resulting veloc- 
ity dispersion (see Heitsch et al. 2009 , for large scale velocity 
di spersion maps in coll iding flows). 

Kle ssen et al.| ([2005) argue that coherent cores can also be 
formed b y gravo-tu rbulent fragmentation of molecular cloud 
material. Kle ssen et al.| ( |2005| ) show the velocity dispersion 
map for some cores m all three projections, and from these 
figures an abrupt increase in velocity dispersion can be identi- 
fied (somewhat in agreement with our observations). How- 
ever, there are important discrepancies between the results 
from Klessen et al. (2005) and the observations: 

1. The increase in velocity dispersion observed in 
NH3(1,1) surrounds the entire coherent dense core, with 
broader lines syste matically found outside the coherent 
dense core. While Klesse rTet al.| ( |2005| ) finds an in- 
crease in the velocity dispersion in more confined re- 
gions (such as a ring around or a stripe next to the co- 
herent core) and with narrow velocity dispersions found 
past these features. 



Fos ter et aT] ( |2009| ) s hows that 81 out of the 83 cores in 
Perseus observed by jRosolowsky et al.| ( |2008| ) display 
subsonic non-thermal motions at their center, while in 
Klessen et aT] ( [2005] ) only a 12-52% of the identi- 



ing mechanism) display coherent subsonic non-thermal 
motions. 

Moreover, it is not clear if any of the models discussed 
above can predict the transition to coherence at densities high 
enough to be observed in NH 3 (1,1). In the case of cores 
formed from shocks this constraint could be extremely impor- 
tant, because the density enhancement generated by the shock 
front can be large enough (a factor of « M 2 , where M is the 
Mach number) to make the detection of NH3(1,1) outside the 
coherent core difficult for highly supersonic turbulence. 

Previous attempts to constrain numerical simulations of 
dense cores using single-pointing survey s of dense gas (e.g., 
Kirk et al.|2007t|Rosolowsky et al.|2008) result in loose con - 



tied objects (which depends on the nature of the driv- 



st raints on simulation s QOlrrier et al.|2008[ |Kirk et al.|2 009). 
In Offner et al. ( 2008]), they present velocity dispersion maps 
derived from synthetic NH 3 (1,1) observations for some cores 
which do not show a velocity dispersion increase similar to the 
one presented in this letter. Clearly these new observations 
allow us to place a different set of constraints on numerical 
simulations that might help to improve the initial conditions 
assumed for star formation. Therefore, it is now the turn of 
simulators to produce synthetic observations from their simu- 
lations that can be compared with those presented here. 
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